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(57) Abstract: A computer-implemented method is provided for valuing and hedging payoffs that are determined by an underlying 
non-marketed variable that moves randomly [200]. The value assigned is that which is obtained by projecting the instantaneous return 
of the future payoff onto the span of marketed assets. An explicit method is provided for determining this value by determining a 
suitable market representative [210]. In a continuous- time embodiment, the methodology is based on an extended Black-Scholes 
equation [220] that accounts for the correlation between the underlying non-tradable asset and marketed assets. Once this extended 
equation is solved [210], the value of the payoff, the optimal hedging strategy [240], and the residual risk of the optimal hedge 
can be determined [250]. In alternate embodiments, the same value is determined as the discounted expected value of the payoff, 
using risk-neutral probabilities for the non-marketed variable. These risk-neutral probabilities are again determined by the relation 
of the underlying variable to the payoff of a most-correlated marketed asset The risk-neutral version of the method applies in both 
continuous-time and discrete-time frameworks, providing asset valuation, optimal hedging, and evaluation of the minimum residual 
risk after hedging. 



WO 2004/008274 



PCT/US2003/021585 



PATENT APPLICATION 

for 

Methods for Assigning a Price to an Asset that is a Derivative 
of a Non-Marketed Variable 

FIELD OF THE INVENTION 

This invention relates to computer-implemented methods within the field of finance. 
More specifically, the invention relates to methods directed at valuing complex un- 
certain cash flow streams and detennining the best way to hedge exposure to such 
streams. 

BACKGROUND 

Frequently it is desirable to assign a price to an asset that is not currently traded, 
or to independently verify a price of an asset that is traded. There are two standard 
methods for assigning such prices (both of which earned Nobel prizes). The first is 
the capital asset pricing model (CAPM) [1] which is a method for assigning a price 
to an asset that has a duration of one period of time (such as one year) and is not 
traded at intermediate times. The second is the Biack-Scholes method [2] for pricing 
financial derivatives, which assigns a price to a derivative (such as an option) on a 
underlying asset that is traded. 

A derivative security is a security whose payoff is determined by the outcome 
of another underlying security. For example, a stock option on a traded stock is a 
derivative, since the final payoff of the option is completely determined by the value 
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of the stock at the terminal time. The Black-Scholes equation is the standard method 
for detennining the value of derivatives. It is based on the fact that, in continuous 
time, it is possible to replicate the payoff of a derivative by a portfolio consisting 
of the underlying security and a risk free asset (such as a U. S. Treasury bill) with 
a fixed risk free interest rate r. The fractions of the portfolio devoted to each of 
its components is adjusted continuously so that the portfolio's response to changes 
in the underlying security will perfectly mirror the response of the derivative. This 
adjustment process is termed a replication strategy. It is argued that the value of the 
derivative is equal to the cost of replication strategy. The cost V is deterrnined by the 
Black-Scholes equation, which gives the cost (or value) V(x, t) for values of x > 0 and 
0 < t < T where T is the terrninal time of the derivative. Specifically, the equation is 

rV (*, t) - V t {x, t) + V x (x, t)rx + \VU*> t)x*o*. (1) 

Here x denotes the value of the underlying security, r is the annual risk free inter- 
est rate, and a 3 is the annual volatility of the underlying security. The notations 
VJ, V X1 14* denote, respectively, the first partial derivative with respect to t t and to X y 
and the second partial derivative with respect to x. The equation is solved with the 
boundary condition V(x,T) = F(x(T))> where F denotes the payoff at time T of the 
derivative. For example, if the derivative is a call option with strike price K t then 
F(x{T)) « max {x(T) -K t 0). 

A more general situation is where a payoff depends on a variable x e but this 
variable is not traded. For example, an option depending on a firm's revenue is of 
this form, because the revenue (which is x c in this case) is not traded. Hence the payoff 
depends on a non-traded underlying variable. In these situations it is impossible to 
form a replicating strategy using the conventional Black-Scholes equation because 
it is impossible to trade the underlying variable. The assumption underlying the 
conventional Black-Scholes equation breaks down. 

Such problems have been studied by other researchers. The idea of using a market 
hedging strategy to rninimize the expected squared error between the final value of 
the hedge and the actual payoff was proposed by Follmer and Sondermann [3], who 
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showed that it was possible in principle. Because minimi'sing the expected squared 
error is equivalent to orthogonal projection of the payoff onto the space of marketed 
payoffs (under a standard definition of projection), the method is often referred to as 
projection pricing. Their analysis, however, is purely abstract and does not exhibit 
any practical method for explicitly finding the hedge by solving a partial differential 
equation or a discrete version of it. 

Merton (in his Nobel Prize acceptance speech [6]) emphasized the importance of 
the problem. He proposed a procedure based on the original Black-Scholes equation, 
but it is essentially an ad hoc method that does not coincide with projection. He 
and Pearson [7] studied a general framework of incomplete markets and formulated 
prices based on a consumer maximization problem. Again their method is abstract 
and does not provide a direct formula for the value that could be used in practice. 

Some practitioners have proposed various ad hoc translations of the standard 
Black-Scholes equation to specific situations. For example, it is common practice to 
artificially increase the volatility of the underlying variable in an attempt to recog- 
nize that the variable is not really a traded asset. These methods are not based on 
optimality, do not fundamentally revise the original Black-Scholes equation, nor have 
any other real theoretical basis. 

An approach that directly addresses the problem presents a solution in terms of a 
market price of risk, which is applicable to all derivatives of a non-marketed variable. 
The market price is difficult to measure, but it has sometimes been estimated [8]. 
However, even if the market price of risk is known, it does not lead to a hedging 
strategy. 

Overall, there has not been an effective and practical method proposed that prices, 
optimally hedges, and computes the residual risk (after hedging) of derivatives of non- 
marketed variables. 
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SUMMARY OF THE INVENTION 

Embodiments of the invention provide computer-implemented methods, based on 
projection, for detennining the value of a derivative of a non-traded variable. In 
addition, some embodiments provide a method for determining the optimal hedging 
strategy and its residual risk. 

Several embodiments of the invention are centered on a new, extended Black- 
Scholes equation of the form 

rV(x e , t) - V t (x ei t) + V^(x Bt t)x& B - - r)] + |^ e (ar e , t^oj, (2) 

where 

Pern = Pem&e/Vm' 

The underlying variable x e is not marketed, (or equivalently, not traded). By this 
we mean, throughout, that it is not marketed as a security that can be held without 
cost and pays no dividend while held. Thus, although oil can be bought and sold, it is 
not traded in this sense because there are storage costs. However, a futures contract 
on oil is traded as a security. In some cases, a variable that is not traded in our 
sense can be converted to a traded asset; (for example, a dividend-paying stock can 
be so converted by subtracting the present value of its dividend). Such cases do not 
require the methods of this invention (although the methods of this invention may 
lead more quickly to the correct result). If our methods are applied to a case where 
the underlying variable is in fact traded, the extended equation will reduce to the 
standard Black-Scholes equation. 

The extended equation (2) uses, in addition to the underlying variable x c , an asset 
x m that serves as a market representative. This representative is itself a marketed 
variable, or a linear combination of marketed assets, in the sense described above. 
This variable enters the equation indirectly through specification of certain param- 
eters of the equation. The constants fjt e and [im are the drift rates of x e and a? m , 
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respectively. p em is the correlation coefficient between x 0 and x m . o% and are the 
annual variances of x e and x mi respectively. 

The replacement of r by - Am (/*m-*0] ^ a coefficient of V ma {x^t)x e in the par- 
tial differential equation distinguishes this extended equation from the conventional 
Black-Scholes equation. This coefficient accounts for the fact that perfect replication 
is not possible in all cases. 

The market representative x m can be chosen in several ways, each of which leads 
to an identical result. In some embodiments x m is chosen to be the Markowitz (or 
market) portfolio of risky assets. According to other embodiments, x^ is chosen to be 
a portfolio of marketed securities most correlated (related) to the underlying variable 
x e . This method is simpler than choosing x m to be the Markowitz (market) portfolio 
and furthermore leads to an optimal mean-square hedging strategy. 

Embodiments of the invention include an implementation based on a solution to 
a continuous-time version of the extended Black-Scholes equation. Thus, according 
to an embodiment, a method is provided for pricing a financial derivative of a non- 
marketed variable represented by x e . The method includes determining a market 
representative that is useful in determining the value of the financial derivative. 
Information associated with the non-marketed variable and the market representative 
is retrieved. A solution to an equation involving a variable V(x ei t) representing a price 
of the financial derivative is then calculated. The equation, which in this embodiment 
is preferably the extended Black-Scholes equation, includes a coefficient involving 
the information associated with the market representative x m and the non-marketed 
variable x c . An output including a calculated price of the financial derivative is 
then generated. Another embodiment includes a solution using formulas based on 
risk-neutral processes. 

The price function can also be found in a discrete-time implementation. The 
underlying process can be discretized directly, or a finite-state model can be used. 
For example, in an embodiment of the invention, a method is provided for pricing 
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a financial derivative of a non-marketed variable represented by x e . The method in- 
cludes retrieving information associated with a suitable market representative x m . A 
solution is then found to a system of equations involving a variable V{x c> k) represent- 
ing a price of the financial derivative at discrete times indexed by k. The equations 
include a coefficient involving the information associated with the variable x m asso- 
ciated with the market representative. An output including a calculated price of the 
financial derivative is then generated. 

In another embodiment, a computer-implemented method is provided for pricing 
a financial derivative of a non-marketed variable represented by a finite-state variable 
B. The method includes detemining a market representative, represented by a vari- 
able A. Risk-neutral probabilities are then calculated using a binomial lattice model 
associated with the non-marketed variable and the market representative. Values of 
a variable V on the lattice corresponding to the variable B are also calculated. The 
calculated values of V represent a calculated price of the financial derivative. 

Embodiments of the invention may also involve computing an optimal hedge and 
other variables of interest. The above embodiments may be easily implemented in 
various different ways on any type of computer, and may be realized as instructions 
stored on a computer-readable medium. 

BRIEF DESCRIPTION OF THE DRAWINGS 

PIG. 1 presents a schematic related to a continuous-time implementation of the in- 
vention. It shows the key variables and how they combine to produce the value of 
the derivative of the underlying variable x e . 

FIG. 2 is a flow diagram showing the steps of the method in the continuous-time case. 
FIG. 3 is a schematic diagram illustrating the projection of V 4- dV onto the space 
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FIG. 4 shows a single step of a binomial model of variable A with probabilities p A 
andl— p^. 

FIG. 5 shows a single step of a binomial model of a variable B when converted to 
risk-neutral form with risk-neutral probabilities of q B and 1 — q B . 

FIG. 6 presents a schematic of the method for the discrete-time case. It shows the 
key variables and how they combine to produce the value of the derivative of the 
underling variable B. 

FIG. 7 is a flow diagram showing the steps of the method in the discrete-time case. 

DETAILED DESCRIPTION: OVERVIEW 

Extended Black-Scholes equation 

Several embodiments of the invention are based on solving an extended Black-Scholes 
equation. This novel equation differs from the standard Black-Scholes equation in 
that it involves a new coefficient in one of its terms. This coefficient allows the 
applications of the equation to be extended to include situations that the standard 
Black-Scholes equation fails to cover. 

A schematic depiction illustrating aspects of an embodiment of the invention is 
presented in FIG. 1. The figure shows an underlying variable x 6i which is specified 
at 100. The variable x e defines a new asset payoff through a function F. Also 
shown in the figure is a market representative x m , determined at 110. The market 
representative x m is extracted from marketed securities, which are determined at 120. 
If this market representative is most correlated with x ei it will depend on z c ; and 
this is indicated in the figure by the dashed line from x e to the market representation 
choice process, 110. The properties of this representative and x e define the particular 
coefficients of a partial differential equation, 130. The function F defines the terminal 
boundary condition at 140. The solution of the equation then yields the value of the 
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new asset at 150. 

Although the new method may at first appear arbitrary since it is based on pro- 
jection, it has a compelling justification. The Black-Scholes approach uses replication 
to render the new payoff redundant, in the sense that the payoff is already embedded 
in the securities market. In a similar way, the new approach approximates the payoff 
and renders the payoff irrelevant in the sense that no risk-averse investor will want to 
own the new payoff because it is inferior to things already embedded in the securities 
market. 

Optimal hedge 

Once the value function V{x Cj t) is found, it is possible to find the portfolio strategy 
that best approximates the new asset, and can be used as a hedging strategy. 

When the hedging portfolio has value H (x ej t), an amount (H(x ei t)—</>) is invested 
in the risk free asset and an amount <j> is invested in the most-correlated marketed 
asset where 

4>{x^t) = V m {x e ,t)x^. (3) 

The minimum variance of the error V(x e , T) — H(x e , T) can be found by solving 
another partial differential equation, namely, 

$(*e,*)+^(*e,t^ = 0, (4) 

with boundary condition S(x ct T) - 0. The value S(x e) 0) is the variance of the 
replication error at T, as seen at t == 0. 

Risk-Neutral Form 

An embodiment of this invention is an implementation as a practical risk-neutral 
valuation method. In this version the drift of the underlying variable x a is changed 

from fi e tO [/i c - Penile - r)]. 



LUB-101/PCT 



8 



WO 2004/008274 



PCT/US2003/021585 



The general formula for value is then 

V(x e ,t) = e- r V-*&[F{x e (T))}, (5) 

where denotes expectation with respect to the risk-neutral model. Hence, in this 
version, the value is simply the discounted expected value of the final payoff, with 
expectation taken with respect to the risk-neutral model. 

The advantage of this form of the valuation is that value can be estimated by 
simulation. Basically, thousands of simulation runs of the underlying variable are 
made and the average discounted value is used as the expected discounted value. 
Such simulation methods are standard practice. A new feature of this embodiment 
is specification of the risk-neutral equation. 

Discrete-Time Version 

The risk-neutral form can be adapted to a discrete-time framework where the dynam- 
ics of random variables are represented by lattices. In the basic model, the underlying 
variable B is assumed at each step to either go up (with probability p B ) or down (with 
probability 1 — pa). A payoff that is function of this underlying variable can be eval- 
uated by using the risk-neutral probabilities, denoted q B and 1 - q B . These, in turn, 
are found by using a market representative as in the continuous-time case. If this 
market representative is denoted A with price (value) v^, the risk-neutral probability 
QB is 

q B = Pb-cov[A^(U b )]\A-v a R]/o* (6) 
= PB-Pitv B)f A&-v A R] (7) 

where 1(U B ) is a payoff of 1 if B is up, and 0 if B is down. The expression "cov w 
stands for covariance, and A(t/ S M 13 defined as cov(l(J/ B ),>4)/variance(i4). 
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A derivative of B is a payoff with a value of G u if up occurs and payoff Q d if down 
occurs. The current price of such a (future) payoff is 

VB-jitofl^+a-tojoj. (8) 

This is the basic method for the discrete-time case. The unique feature is the 
formula (6) for the risk-neutral probability. 

Recursive Solution and Hedging 

If the payoff occurs at the end of some finite number of steps (at the end of a lattice 
structure for B), the value function can be found recursively. The recursive solution 
is 

= ifaaVkfoU) + (1 - te)V5,(4-i)]. (9) 
where denotes the upper successor state to s k -x and sjj_ x denotes the lower 
successor state to s fc _i. The process is started with the tenninal boundary condition 
specifying the payoff of the derivative G. 

An optimal hedge H can be found. If at the beginning of time period A; the hedge 
portfolio's value is H^\ t then H k _ x - yV k ^ is invested in the risk free asset and 
7Vt-i is invested in the most-correlated asset A, where 

7 = cov(V^-i>^)/"i (10) 

The variance associated with optimal hedging can also be found by a recursion, as 
discussed later. 

A Complete Methodology 

Embodiments of the invention provide complete methodologies for treating payoffs 
that are derivative of non-traded variables. It has a strong justification and a straight- 
forward representation. For many problems, the new method can be put into practice 
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by augmenting standard methods of financial computation. The requisite equations 
or lattice structures are similar to those derived from the Black-Scholes methodology 
and hence the computational methods are similar. These are methods for solving 
partial differential equations, methods for evaluating lattices, and simulation meth- 
ods. A significant additional step of the new method is determination of the market 
representative x m , and estimation of the constants /* e , ^ (the drifts), the variance 
of the market representative and the value of ftm- 

Applications 

Embodiments of this invention have applications in numerous financial areas. For 
example, the methods may serve four functions. The first is pricing, as in detenniri- 
ing the fair value of an option on a non-marketed variable. The second is design, as 
in the design of optimal contracts, project plans, business arrangements, and various 
agreements. The third is hedging, to minimize the risk associated with a non-traded 
asset by use of offsetting market participation. The fourth function is that of risk 
assessment, determining the residual risk after optimal hedging. In all cases, a finan- 
cial payoff is tied to the performance of a variable that moves randomly but which is 
not a marketed financial security. The variable that defines the payoff is termed the 
underlying variable. 

The key idea is that when the underlying variable of an asset is not marketed, and 
hence the standard Black-Scholes equation does not apply, a market representative 
can be used instead. 

One large area of application is to the pricing and hedging of "off-exchange" 
derivatives. For example, it may be desirable to hedge grapefruit production. Here the 
underlying variable is grapefruit price but there may be no financial instrument (such 
as a grapefruit futures contract) that is directly related to this underlying variable. 
Hence grapefruit production is not a derivative security and the Black-Scholes theory 
is not applicable. Nevertheless, a combination of orange juice futures contracts may be 
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most-correlated with grapefruit prices and thus can serve as the market representative. 
This representative can be used (according to the valuation equation) to find value 
and also serve as the basis for an optimal (but imperfect) hedging strategy. Similar 
situations arise with off-exchange trading of energy contracts, agricultural products, 
metals, foreign exchange, and many others. 

One important example is that of real estate projects. Typically, there is no 
underlying security for a real estate project, and thus the Black-Scholes equation 
does not apply. However, the payoff of such a project is likely correlated with interest 
rates, real estate investment trusts, and regional economic variables. A combination 
of these can serve as the appropriate market representative and hedging strategy. 

Another general area of application is the determination of long-term futures 
prices. Typically, futures contracts extend several months, or at best a few years. Yet, 
many large projects span several years and correct pricing and hedging is important. 
The methods of this invention can be used to obtain such prices and hedges. 

Bonds are often issued that have unique risk, not completely correlated to mar- 
keted securities. For example, a corporate bond's price may be only loosely correlated 
to interest rates and to the stock price of the associated firm and to the stock prices 
of firms in the same or related industry. Hence the bonds are not perfect deriva- 
tives, but the methods of this invention can be used to evaluate them. Large physical 
projects, such as oil rigs, dams, various other infrastructure projects and private or 
public works have cash flows that are to some extent associated with risks that are not 
derivatives of market securities. Projects in foreign countries have associated country 
risk, which can only partially be hedged with bonds and insurance. 

Contracting is another large area of application. For example, an electronics firm 
may plan to purchase a great deal of DRAM during the year. The firm can design 
contracts with options and price caps that provide reduced risk at favorable prices. 
The design of such a contract can be deduced from the market representative and 
the methods of this document. Contracts for movie and publishing rights can also be 
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designed this way. 

The purchase of intellectual property rights through licenses is a fertile area of 
application. For example, a pharmaceutical firm may wish to license a biochemistry 
patent from a university, and both the design and the pricing of such a license could 
be carried out with the methods. In this case, the market representative may be a 
combination of the security prices of other pharmaceutical companies, HMOs, and 
other insurance companies. 

Insurance and various guarantees generally embody a great deal of risk. Such risk 
can be priced and hedged by the insurance company or guarantor by these methods. 
The risk of fire in an area may, for example, be correlated with weather derivatives 
and energy prices. 

Mortgages contain interest rate risk, but also prepayment risk. The prepayment 
risk may be correlated (imperfectly) to stock market indices, real estate prices, as 
well as interest rates. 

Real options is a general term which includes options related to many of the exam- 
ples listed above, but especially to business projects. The methods of this document 
can improve both the evaluation and design of such projects by using a theoretically 
justified valuation formula rather than an ad hoc use of the Black-Scholes equation. 

Another important area of application is that of overall risk assessment. The 
residual risk of each individual project or derivative can be determined by the variance 
given as the solution to a partial differential equation that is used in embodiments of 
the invention. The risk associated with a collection of such projects can be found by 
solving a similar equation that depends on the value function for each of the projects 
separately. 



LUB-101/PCT 



13 



WO 2004/008274 



PCT/US2003/021585 



DETAILED DESCRIPTION: 



CONTINUOUS-TIME CASE 



This section presents example application areas and gives a detailed description in 
both the continuous-time and discrete-time frameworks. 

Outline of Steps: Continuous-time framework 

An outline of the steps used in one embodiment of the invention are indicated in FIG. 
2. 

1. Set Up. In this step, 200, we identity an asset of interest which gives a future 
payoff (or payoffs) according to the value of some variable that moves randomly 
with time. This variable is termed the underlying variable. If this variable is 
itself traded (and is governed by geometric Brownian motion (GBM)), then the 
asset of interest is a pure derivative and the standard Black-Schoies equation 
can be used. If the underlying variable is not traded, then the extended Black- 
Scholes equation is used. 

For example, an asset might be the future harvest of a grapefruit farm. The 
underlying variable determining payoff is grapefruit price which varies randomly 
from the current time until harvest. This underlying variable cannot be traded. 
It is true, of course, that grapefruit can be bought and sold at any time, but 
grapefruit cannot be held as a security since there are storage costs and grape- 
fruit is perishable. The grapefruit price is denoted x e and is assumed to follow 
a GBM process of the form 



where z c is a standardized Wiener process. The process starts at time 0. The 
asset to be priced has a payoff at time T of F(x e ,T). 



dx e = fi e x e dt + a e x e dz B 



(11) 
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There axe n (perhaps thousands of) securities that can be traded. Each of these 
follows a similar processes, namely, 

dxi s= fiiXidt + (XiXidzi (12) 

where x< is the value of the i-th asset, for 1 = 1,2,..., n. The standard- 
ized Wiener processes Zi (of zero mean and unit variance) axe correlated, with 
cov(dz*, dq) = pijdt Among these n assets there is an asset that is risk free at 
each instant. By forming an appropriate combination of the marketed assets, it 
may be assumed without loss of generality that the risk free asset is x n . Then 

dxn = rx n dt. 

In the standard Black-Scholes method, only the underlying security, the risk free 
rate r, and the payoff function is specified. No consideration is given to other 
securities, because the underlying security is all that is needed. In embodiments 
of this invention, other securities are considered in order to select a market 
representative, as explained in the next step. 

2. Deterniination of a Market Representative. The extended Black-Scholes equa- 
tion uses a market representative x mj which is a combination of marketed assets. 
In this step, 210, this market representative is determined. Generally, this mar- 
ket representative can be chosen in one of several ways. 

(a) A most-correlated asset is a marketed asset x c (or a combination of mar- 
keted assets) whose instantaneous return is most correlated (or close to 
most correlated) with dx G . In theory, it is obtained by solving the equa^ 
tions 

n-l 

X>« a i = *fa (13) 

for each i = 1,2, . . .n — 1. The resulting o^'s are then normalized to sum 
to 1. The most-correlated asset is the asset x c with instantaneous return 

dx„ 2=J dxi 
a*c i=\ Xi 
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x c can be used as the market representative x m . 

In practice, a most-correlated asset is obtained by considering only those 
securities that axe obviously closely related to the underlying variable x c . 

In the case of grapefruit, for example, a most-correlated asset might be a 
futures contract on orange juice. (There axe no grapefruit futures.) Such 
contracts are marketed securities, and are likely to be highly correlated 
with grapefruit prices. 

As another example, in seeking a most-correlated asset for DRAM (com- 
puter memory chips), one might consider the stock of several companies 
that produce DRAM and then find the combination of these stocks that is 
most correlated with DRAM prices. 

Ideally, this step entails identification of suitable candidates, evaluation of 
the covariances of these candidates among themselves and with the under- 
lying variable, and optimization of the combination of these candidates to 
provide the highest correlation. These steps are not difficult provided one 
has historical data on the security prices and the underlying variable. 
However, it should be understood that correlations, variances, and drift 
rates are measured approximately, and hence in practice exact parameters 
are not known, even for the standard Black-Scholes equation. Likewise, 
it is impossible to prove that a given asset is, in fact, the asset that is 
'most 1 correlated with the underlying variable. As in any method driven 
by estimates from limited data, one attempts to balance accuracy with 
data gathering and processing effort. Thus, in the context of the present 
description, a 'most correlated' asset or variable is defined to include not 
only the maximally correlated asset or variable, but also approximations 
thereto. 

Although in practice a 'most-correlated' asset is not known with certainty, 
this embodiment is b^tinguished from the following one by the fact that 
the search for the strongest correlation is restricted to obvious candidates 
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within one or a few related market sectors or industries. For example, the 
steel industry likely would not be considered as a candidate for an asset of 
grapefruit production. This contrasts with the Markowitz portfolio (gen- 
erally approximated by the market portfolio) which includes all marketed 
securities. 

(b) The Markowitz or market portfolio. This alternative embodiment produces 
the same price as the method described above. This method is typically 
less convenient than the first and does not lead to an optimal hedging (or, 
equivalently, replicating) portfolio. In theory, this alternative is determined 
by solving the system of equations 

These a^'s are then normalized to sum to 1 and used to define the portfolio 
with instantaneous return 

dx M S^i dx, 

Then x M can be used as the market representative x m . 
As an approximation to the Markowitz portfolio, the market portfolio con- 
sisting of the capitalized weighted average of the marketed assets may be 
used. For example, a diversified mutual fund can be used, or the S&P 500 
index. Each of these is nevertheless only an approximation. It would be 
essentially impossible to determine the true Markowitz portfolio because 
of the inherent error associated with deterrmning the required parameters. 
In contrast to the previous embodiment, the steel sector, and indeed all 
sectors of the economy, would be represented in this method. The ad- 
vantage of this method is that the Markowitz (or market) portfolio, once 
found, can be used for derivatives of any underlying variable. 

(c) Local index. Suppose there are two non-marketed variables x e i and X& 
that determine various derivative assets. Let x^ and X& be corresponding 
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most-correlated risky marketed assets. Next let x m be the Markowitz com- 
bination of Xd and x<&. Then the single asset x m can be used as a market 
asset for all derivatives of x cl and x^. This procedure can be directly 
extended to any number of underlying variables. 

To complete the set up, the payoff function F(x £ (T)) ) the parameters p c and 
a e of the underlying variable, the parameters /i™, ft™, and the interest rate 
r are specified. As stated earlier, in practice it is not necessary to specify 
the parameters of every marketed security. A market representative can be 
determined by consideration of those assets closely related to the underlying 
variable or, alternatively, to a broad index of the market. 

3. Define the Differential Equation. In step 220 various terms are specified for the 
extended Black-Scholes equation: 

rVfo.O ( 15 ) 
where 

Pem = EfdayfenJ 

The subscript m denotes the market representative. The boundary conditions 
are 

V(x e ,T) - F(x e (T)) 
y(o,t) = e-^-'^o). 

Equation (15) plays a central role in several embodiments of the invention. It 
differs from the standard Black-Scholes equation in the coefficient [^ e — Pem[t*m— 
r)]. This coefficient replaces r in the standard Black-Scholes equation as a 
coefficient of T£ 0 (x e ,t)x e . The new coefficient uses properties of the market 
representative x^. Equation (15) reduces to the Black-Scholes equation if the 
new asset is a pure derivative so that x m = x* is the most-correlated market 
representative. 
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4. Solve the Differential Equation. In step 230 the differential equation is solved. 
This can be done analytically in some cases. Most often it will be solved by 
standard numerical techniques. Once the parameters are known, solution of 
the equation is no more difficult than solution of the standard Black-Scholes 
equation. It is only one coefficient that is different. 

The value of V(x e (0), 0) is the proper price of the derivative. The procedure can 
be terminated here, with the proper price. Additional information is obtained 
in the next two steps. 

5. Find the Optimal Replicating Portfolio. Optionally, step 240 calculates the best 
replicating portfolio. One may replicate the market portion of the risk in the 
asset by holding an appropriate portfolio of market securities with initial value 
equal to V(x e (0), 0). This portfolio is a combination of the most-correlated asset 
and the risk free asset. The amount to be held in the most-correlated asset as 
a function of x e and t is 

where in this case, we use the subscript c to emphasize that the market rep- 
resentative is taken to be the most-correlated asset. Holding the negative of 
this replicating portfolio produces an optimal hedge for the derivative. This 
procedure is essentially identical to that used in standard derivatives of the 
Black-Scholes type. The difference is that in the Black-Scholes case the hedge 
is perfect, while in embodiments of this invention the hedge may not be perfect. 

For example, in the case of grapefruit production, the optimal hedge would 
be a portfolio consisting of the risk free asset and a position in orange juice 
futures. The proportions of the two components of the portfolio would be 
adjusted frequently according to the formula for ^. If the hedge is in place, the 
hedge portfolio yields cash immediately equal to V(x e (0) t 0) and the final net 
payoff (consisting of sale of the grapefruit production and settling of the hedge 
contract) is on average zero, but it has a residual variance. 

Ideally, the portfolio adjustment would be continuous, but in practice periodic 



LUB-101/PCT 



19 



WO 2004/008274 PCT/US2003/021585 



adjustments (daily, weekly, or monthly) are sufficiently accurate. These types 
of adjustments are common in pure derivatives methodology as well, where the 
ideal hedge is perfect. 

6. Find the Error Variance. Optionally, step 250 calculates the tracking error. 
The variance of the part of the payoff that cannot be hedged is found by solving 
another partial differential equation which uses the solution to (15) as an input. 
The equation is 

t) + S^x e + 4- e*< r -<>[K.a e s c ] 2 (l ~pl) = 0 (16) 

with boundary condition S(x e9 T) = 0. The value of S(x e (0), 0) is the variance 
of the error at time T as seen at time 0. This may be extended to computation 
of the variance associated with combinations of derivatives. 

For example, in the case of grapefruit, the value of 5(x c (0), 0) gives the variance 
of the exposure at time T with an optimal orange juice futures hedge. That 
is, with the hedge, V(x e (ti) t 0) is obtained at t = 0 and the net amount (of 
the combined grapefruit production and hedge) attained at t = T has expected 
value 0 but variance equal to S(x e (0), 0). 

This calculation of error variance has no analogy in the standard Black-Scholes 
case because in that case the error is always zero (in theory). Hence, this 
method for finding the hedge variance is new and important. It is not difficult 
to solve the appropriate partial differential equation, for indeed, it is similar in 
structure to the equation for value. 

Extensions 

There are several important extensions of the method. For example, 

1. Inclusion of cash flows that occur with time, with incremental cash at time t 
being of the form h(x ey t)dt. 
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2. Varying parameters. The parameters of the model can vary with i, x e , and the 
Xi's. In general this leads to a partial differential equation of higher order. 

3. Parameters, such as the risk tree rate, may be governed by stochastic processes. 

4. Additional variables. Additional non-marketed variables can be introduced. 
These serve as "state variables" for the system. For example, a state variable 
might be total industry productive capacity, or estimates on the probability 
that certain legislation will be passed. 

5. Estimation variables. A suitable non-marketed but observed variable may be 
the estimate of an unobserved variable that serves as the underlying variable for 
the payoff. If the estimate converges to the actual value at the time of payoff, 
the estimate may be used at all points instead of the original variable. For 
example, the best estimate for yearly revenue may converge to the actual figure 
as the year ends. 

6. Non-market random components. If the final payoff is a function of the form 
F*(x et y, T) where y is random and independent of x e and the market, then we 
define F(x et t) = Ey[F*(x CJ y,T)] where Ey denotes expectation with respect 
to y. The value of the payoff can then be found with this F and it will be 
the proper price. In this case, the boundary condition for auxiliary differential 
equation for variance is 5(x e ,T) = variance y [F*(rr c , y,T)]. 

7. Alternative Processes. The method of relating the non-market variable to the 
market in order to obtain the correct price can be extended to alternative pro- 
cesses, including jump processes, and moving average processes. 

8. Market Cash Flows. In some cases the cash flows or payoffs may depend on 
marketed variables as well as non-marketed variables. In such a case the value 
function V will depend on both x e and the marketed variables. The most- 
correlated asset for the marketed variables will be those variables themselves, 
while the market representative for the non-marketed variable is found as in the 
basic case. 
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9. Path Dependent Cash Flows. These are cash flows such as max^O, x c (t)], 0 < 
t < T that depend on the actual path taken by x e rather than on its instan- 
taneous or final value. The risk-neutral process can be used to evaluate such 
situations by taking the risk-neutral discounted expected value of all cash flows. 

Brief Derivation of the Equations 

We refer to step #1 (Set up) to define the problem. Hence, we have an underlying 
variable x c governed by geometric Brownian motion (GBM) over 0 < t < T as 

dx e = (i c x e dt + a c x e dz C} (17) 

and likewise there are n marketed assets that also are governed by GBM. There is a 
special asset of interest whose payoff at time T is F(x c (T)). 

We propose an (initially unknown) price function V(x ey t) on the time interval 
[0,jT] with terminal value V(x Ci T) = F(x e ). It will be defined by instantaneous 
discounted projection. 

At a fixed time t we define the space M as the linear space generated by all instan- 
taneous marketed returns (dx*)/^. That is, M is made up of all linear combinations 
of those instantaneous returns. Symbolically, the value of our V(x e) t) will satisfy 

V(x ei t) = P{V(x cy t) + dV(x ei t)|M}. (18) 

which denotes the discounted projection P onto the space M of instantaneous mar- 
keted returns. Specifically, V(x ej t) is found by: first calculating the payoff m* at 
time t + dt in M that is closest to V(x e> t) + dV (x c , t) in the sense of minimizing the 
expected square of their difference; and second, computing the price at time * of m* 
(which is defined by linearity since it is a combination of marketed returns). Briefly, 
we say that the value at t of the payoff at t + dt is given by discounted instantaneous 
projection. This process of projection onto M is illustrated schematically in FIG. 3. 

Since payoff V(x et t) at t -f dt is certain, its price at t is the discounted value 
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(1 - rdt)V(x e) t). Thus (18) becomes 

rV(x ej t)dt = P{dV (x e) i)|M}. (19) 

We shall only keep terms that are first-order in dt. Hence, in particular, P{di} = dt. 
Substituting the Ito formula 

dV(x Cl t) = [V5(x e> t) + V^(x ey t)^x c + %V XeXc oixl]dt + K.<wl* (20) 

into (19) and keeping only first-order terms in di, (19) becomes 

rV (x e , t)dt = [K(x e> t) + ^(Xe, *Kx c + ^.of x*]d* + V^<7 c x c P{d* e |M}. (21) 

The projection of dz t is a linear combination of rdt and the differential of a 
marketed asset most correlated with dx c . A marketed asset x e with instantaneous 
return dx c /x c most correlated with dx e /x e satisfies dx c = [icX e dt + o<x c dz c where 
/ic and a c are each formed as (identical) linear combinations of the /Vs and a^s, 
respectively. z c is a standardized Wiener process. 

Standard methods show that the projection of dz* onto M is Pecdzc. In addition, 
from 

x c = P{x c 4* /icX c di 4- x^dzc} 

it follows that 

P{dz c ) = ^""^ dt. (22) 
Substituting this in (21) produces 

rV(x Ci t)dt = [V t (x e ,t) + V^(x c ,i)M^^ 

= [V t (x et t) + V x .(x C7 t)x e [p e -A.fo.-r)] +iVx e xe(^,t)^]d*, 

where we define 

» 

Canceling the dt we have the final result 

rV(x Bl t) - y t (x c ,i) + ^(x e ,i)x c [ Me -^ cc ( Mc -r)] + |W*,i)aJa£. (24) 
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This is the basic extended Black-Scholes equation. 

Once equation (24) is solved with the boundary conditions V(x ey t) = F[x e ), 
V(Q,t) = e- r ( r_t >jP(0), the value V(x c (0),0) k the value at time 0 of the derivative. 

Market Representatives 

Details of how the two basic alternative market representatives, a most-correlated 
asset and the Markowitz portfolio, are computed are given here. In practice, these 
detailed calculations are greatly simplified (approximated) by use of intuition or equi- 
librium arguments, as mentioned below. 

1. Most-Correlated Asset. 

Arrange the n marketed assets so the first n — 1 are risky and the n-th is risk 
free. The return of a most-correlated asset is of the form 

At nr ~ 1 dx- 

where the a^s do not necessarily sum to 1. 

To maximize the correlation of this return with dx„, the o^'s solve 

_ / dx c dx e \ 

max E 

\x c x e J 



fdx c \ 

subject to var f — J < 1. 



Introducing a (positive) Lagrange multiplier for the constraint, and representing 
dx c /x c in terms of the a^s, this is equivalent to solving 



max 



it— 1 n— 1 

£ eyerie - A <Xi<Tij<Xj 

»=1 ij 



(26) 



Since any solution can be scaled by an arbitrary positive constant and still 
preserve the property of being most correlated, without loss of generality one 
may set A = 1. Then (26) becomes 

£ <W = <r rf (27) 
4=1 
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which is easily solved. 

The current price (at t) of dx c /x c is ££1 To make the price equal to 1, the 
oVs are scaled so that they do sum to 1. Then the correct x c is given by (25) 
with the scaled a/s. 

In practice, a market representative may be chosen as an asset that is obviously 
closely related to the variable x e . For example, an orange juice future is closely 
related to grapefruit price. It is probably not necessary to consider oil and 
metal futures or the thousands of marketed stocks to find a closely related 
market representative. A refined approach is to consider a family of agricultural 
futures and, based on historical returns, compute the most-correlated portfolio 
from this family. The selection of a closely related market representative is, in 
practice, partially art and partially data collection and computation. 

2. Markowitz portfolio. The Markowitz portfolio xm has the advantage that it can 
be used to price a derivative of any process x e , but its disadvantages are that 
is difficult to use, less intuitive, may not exist, and does not lead to an optimal 
replication. (The use of this portfolio in the single-period case is described in 
Luenberger [10].) The Markowitz portfolio x^ is the portfolio of purely risky 
assets (the first n — 1 assets) that has price 1 and maximizes 

If it exists, the set of o»/s that achieves it satisfies 

E^ = cU-r) (28) 

1=1 

for the constant c that makes the Oj's sum to 1. This portfolio is an alternative 
for purposes of valuation only, not for determining optimal replication or hedg- 
ing. In practice, a broad market index, such as the S&P 500, is used as a proxy 
for the Markowitz portfolio. 

i 

3. Local index. Suppose x e i and x c2 are non-marketed underlying variables and 
x& and X& are the respective most-correlated marketed assets. The Markowitz 
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combination of x& and x& is the portfolio x m = ai^d = a%x<a where the weights 
ai and a 2 are determined as in the above discussion of the Markowitz portfolio 
under the assumption that x& and r are the only marketed assets. This x m 
is a market representative for any derivative of x<& and ar c2 . It serves as a local 
index, and may in practice be the index of an industry to which x el and X& are 
related. This local index can be used to value a variety of derivatives. It cannot 
be used as the basis for optimal hedging, but it can be used to form approximate 
hedges. The method extends to any number of underlying variables. 

Important Special Case. 

If the payoff is a put or call option on the variable x et a closed-form expression for the 
solution to (15) applies, which extends the special case solved by Black and Scholes. 
If the strike price of the call option is K t the payoff is max[rr e — K y 0]. The value of 
the call option is 

V(x ei t) = x e e^- r ^ T ^N(d 1 ) - Ke-rV-WNith) (29) 
where N(d) denotes the value of the standard normal distribution at amount d t and 

da = d 1 -a e y/T=i (31) 

U> = fr-pem(frn-r) y (32) 

where is the drift rate of the market representative x m and j5 m is the beta of x c 
and Xm. 

The value of a put with strike price K can be found from an extended putr-call 
parity formula 

C - P + Ke-* = x e e^ T , (33) 
Here C is the price of the call, and P the price of the put. 
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Optimal Replication 

It is possible to optimally replicate the final payoff by trading in the marketed assets. 
The replication is imperfect, but the error has zero expected value and is uncorrelated 
with all marketed assets. This is, therefore, the best that can be done within a mean- 
variance framework. The trading strategy for replication can be derived from the 
value function. 

As in (20), the value function follows the Ito process 

dV(x c ,t) = [V t (x c , t) + V Xa (x e> t)fJL e X e + |K.x.(Xc*)^]dt + V^x.dz.. (34) 
One may write 

dZe = PeedjZc + yjl - fadZp, (35) 

where dzp is a standardized Wiener process uncorrelated with all marketed assets and 
where p w = ^/(a^) — PecPc/<7e- With this substitution and substitution of the 
term ^Vx^ Xo (x e} t)o*x* from the extended Black-Scholes equation (15) it is possible to 
transform (34) to 

dV (x e , t) = [{V(x et t) - <f>)r + # c ]dt + <fxr c dz c + Sdz,, (36) 

where 

<l>(x ei t) = V^x^xJec (37) 

S(x et t) = V^x^o^yjl-fc. (38) 

Motivated by (36), the process H is defined by 

dH(x ai t) = [(H{x et t) - <t>)r + tajdt + tfocdzc, (39) 

with initial condition H(x cy 0) = V(a: e , 0). 

From (39) it is clear that dH(x ei t) is in M at (x e) t) and that, overall, H(x e) t) is 
generated by linear combinations of marketed returns. From 

d{V-H)=r(V-H)dt + 6dz p . (40) 
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and the initial condition V — H = 0, it is clear that at time 0 there holds E[V(x e , t) — 
H(x Ci t)] — 0 for all 0 < t < T. Furthermore, since V — H is a linear integral of 
d^'s, it is orthogonal (that is, uncorrelated) with all marketed returns. Thus the 
random variable H(x e> T) as seen at time 0 is the projection of V(x ey T) onto the 
space of marketed assets over 0 < t < T. 

It follows that the optimal replication is simply H y governed by (39). This repli- 
cation requires only initial cash of V(x Ci 0) and no other infusions or withdrawals. At 
each instant H — (/> is invested in the risk free asset and (j> in the most-correlated asset 
(for a total of H). Holding — H serves to optimally hedge the derivative. 

Projection Error 

The variance of the residual replication error V — H at T can be found as a solution 
to a partial differential equation adjunct to the general pricing equation. 

Define the difference variable D - V - H, and let U = e^'^D. Then from (40) 



where 6*(x t1 t) = e r P-*>6(x c ,t). 

Let S(x ey t) be the variance of U(T) when in (41) U is initiated with U = 0 at the 
point (x e , t). This variance is 



= **(a; e , tOd^, 



(41) 




(42) 



where E denotes expectation at t. 



We write (42) as 




| J dt + E^ t 5*(x e (a),s) 2 ds 
\ 2 dt + E[S(x ei t) + dS{x e ,t)]. 



(43) 
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Hence, 



E[d5(x e ,t)]+<J*(x CJ t) 2 dt = 0. 



(44) 



Expanding (44) using Ito's lemma, we find 



St{x e ,t) + S^x, + \t^M + e*<r-*\y Xo (xe,t)a c x e )*(l - pi) = 0, (45) 



with boundary condition S(x 6i T) = 0. The value S(x ei 0) is the variance of the 
replication error at T, as seen at t = 0. 

Universality 

As a by-product of the optimal replication equations, we obtain a compelling justi- 
fication for methods used in embodiments of this invention. When the new asset is 
priced according to teachings of this invention, the optimal replication H has the same 
expected value as the asset with payoff F(x t (T)). In addition the optimal replication 
has lower variance than the asset since V(x e ,T) = H(x Cl T) + error, where the error 
is uncorrelated with ail marketed assets. This means that every risk-averse investor 
will prefer the replication over the new asset. This property is termed universality, to 
reflect that it is true for everyone. Another way to state this property is that if the 
new asset is priced according to methods of this invention, every risk-averse investor 
will elect to include the asset only at the zero level (preferring neither to be long or 
short the asset). In this sense, the price renders the asset irrelevant in the market, 
for it is dominated by the best replicating asset which is already available. 

Risk-Neutral Process. 

These methods can be transformed to formulas based on a risk-neutral processes. 
The novel feature is that the risk-neutral process is defined in terms of a market 
representative x m defined earlier. The appropriate risk-neutral process is 



cb e = cox^dt -f a c x £ dz e 



(46) 
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where 

" = Me-/M/im-r). (47) 
In these terms the proper price of the payoff F(x e (T)) is 

P = e-*r$[F(x*(T))) (48) 

where E denotes expectation at time 0 with respect to the risk-neutral process. One 
way to solve such a problem is by simulation- Many paths of x e are generated ac- 
cording to the risk-neutral process and the resulting terminal payoffs are averaged to 
obtain an estimate of the risk-neutral expected value. The discounted value of this 
estimate converges to the correct price as the number of simulation runs increases. 

A significant advantage of the risk-neutral process is that it can be used to 
value path-dependent cash flows; where the payoff depends on the particular path 
taken by x e as well as its final value. For example, the Asian option with payoff 
max[maxo<t<rx e (r.) - K t Q] can be evaluated using the risk-neutral process. The 
proper price is again 

which can be evaluated most easily by simulation, but also by special grid methods. 

DETAILED DESCRIPTION OF THE INVENTION: 
DISCRETE-TIME CASE. 

Approximation 

One way to work in a discrete-time framework is to directly discretize the model and 
the partial differential equations. For example, the new version of the process for x e 
becomes 

x e {k + !) = (! + /ieAr> e + <JeX e e(k)y/Ai. (49) 
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In this equation At is the length of the time step of the model and e(fc) is the value 
of a standardized normal random variable uncorrelated with previous or future such 
variables. The market variables are transformed to discrete form in a similar way. 
A corresponding discrete-time version of the extended Black-Scholes equation can be 
developed. It is 

W - 1)) - jg [mWM - cov(v * (x » (fe) - «"(*>>PM*)1 - *w)j (50) 

where R = e tA< and where x m is the market representative (either the Markowitz 
portfolio or a most-correlated market asset or a combination of these.) The expected 
values are taken at time fc-1. The quantity is the price at time k-lofx mi is 
the variance of x m (k) as seen at A: — 1. The formula is valid for any discrete process 
that governs x c , although the market variables follows a fixed process as in (49) (but 
with i's instead of 6 as subscripts). The extensions applicable in the continuous-time 
framework are also applicable in discrete time. 

The discrete-time approximation (49) can be converted to risk-neutral form as 

x e (k + 1) = (14- u>At)x e + <r c x e €(k) VS, (51) 

where u = /x c - f^fa* - r). 

Finite-State Models 

A single stochastic variable A is represented by the two-state model shown in FIG. 
4. The Ua outcome is considered an "up" move, and the D A outcome a "down" 
move. The model has three degrees of freedom (one probability and two node values) 
and hence the expected value and variance of the variable can be matched with one 
remaining degree of freedom. The model can be extended over several time periods 
as a binomial tree or (frequently) a binomial lattice. 

There are two standard approaches to parameter matching: additive and geomet- 
ric, corresponding to matching moments of the variable itself or matching moments 
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of the logarithm of the payoff. Geometric matching is a natural choice for processes 
governed by geometric Brownian motion. Additive matching is frequently used for 
discrete-time models. A binomial model for one of these is easily converted to a cor- 
responding binomial model for the other. For small At defining the time step of the 
model, the two approaches are nearly identical. 

Suppose that A and B are variables as above, and suppose G is a function of (a 
derivative of) B, G is defined by the two values G u and G d corresponding to whether 
Ub or Db occurs. There is an important relation between the covariance of such a 
derivative with A and the covariance of B with A that holds when B is described by 
a binomial model. This is spelled out in the following result, easily proved by algebra: 

cov(AG) = ^^cov(AB). (52) 

Suppose A is the market asset most correlated with B. The above result guaran- 
tees that in the case of binomial models, A is also most-correlated to all derivatives 
of B (since the two covariances are proportional). It follows that the three-variable 
model (with A, B t and the risk free total return J?) can be used to price all derivatives 
of B. 

Paralleling the development for the continuous-time case, A is a market represen- 
tative and it may be taken to be a marketed asset most-correlated with the underlying 
B or, alternatively, as the Markowitz (market) portfolio. 

Let G be a derivative of £, defined by its two values G u and G d . The projection 
price is (with an over-bar on a random variable denoting expected value) 

v G « i[E(G)-cov(G, A)(X-v A R)/o* A ] 

- -g [pbG u + (1 - PB )G d - [cav(l(Efc), A)G U + cov(l(D B ), A)G d ]& - v A R]/* 2 A ] 
= ^ { [pb - / VbmF - v A S\] G u + [(1 - PB ) - Ap fl P - vaR}] G d } , (53) 

where 1(U B ) and 1{Db) denote payoffs of 1 if U B or D By respectively, occurs; and 
where (h& B ) t A = cov[l(tf B ), A]/o\. 
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This can be written as 

Vo = 5 [QbG u + (1 - to )Gf j (54) 

where 

qB=PB- P\{u B )^ P - tu#] . (55) 
Expanding the beta term, we can calculate explicitly 

— Pb Pb{1 - Pb)[Au - 4fp - tU-R]/<7^ (56) 
where 4* = E[A|i7 fl ] ^ = E[A\D B ]. 

This is the basic method for the discrete-time case. The unique feature is the 
formula (55) or (56) for the risk-neutral probability of an up move. Once this proba- 
bility is determined, it may be used in place of the true probability p A for the purpose 
of evaluating payoffs that depend on B. In other words, for purposes of evaluation 
the binomial model takes the form of FIG. 5. 

Recursive Solution 

The single-period structure can be extended to a multiple-period framework by piecing 
together single periods, and this leads to a recursive solution. 

Suppose the nodes of the lattice are numbered by the time index k and the state 
index s k which is the level of the node counting from the bottom. Briefly, we write 
the node as (fc, s k ). The value function of the derivative is a value at each node, and 
described as V k (s k ) at time point k and state s k . 

The recursive solution is 

IW**-!) = gfeBliMU) + (1 - q B )V k (8tx% (57) 

where denotes the upper successor state to s k ^ t and denotes the lower 
successor state to s k -i. The process is started with the terminal boundary condition 
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specifying the payoff of the derivative G. If there are additional payoffs along the 
way, they are incorporated step by step in the usual manner. 

This method can be extended in the same ways as the continuous-time version. 

FIG. 6 is a schematic depiction of the method for the discrete-time case. This 
schematic has the same basic structure as that of FIG. 1 which depicts the schematic 
for the continuous-time case, with the exception that the variable V is found through 
a recursion process 600. The figure shows how the underlying variable B defines the 
new asset and how the market representative A is extracted from the market. The 
properties of these two variables define qs which defines the recursion, leading to the 
value of the new asset. 

Error Process 

Denote by {T4|M} the projection of V k on the market at time k — 1. This will be the 
projection of V k onto the space spanned by R and a market representative A most 
correlated with B. It is easily shown that 



where e k is uncorrected with the market. Using the pricing equation (50), one has 



{V k \M} = V h 4- cov(^, A)(A - 



Hence 



V h = V k + cov(V k> A)(A - + 



(58) 



BV k -i =Vk + cov(V fe , A)(R - 3)/o* . 



Eliminating V kt (58) can be written as 



V h = JH^i - cov(V fc , A){R - 3)M + cov(V fc , A)(A - 3)/cfl + e k 
= [(l-7)* + 7A]VU + £* 



(59) 



where 



7 = cov(V^-i>^)M- 
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In the same manner as in the continuous-time case this shows how to select the best 
approximating (or replicating) portfolio H, with an amount H k _ x - *yV k _ x in the risk 
free asset and 7^ in the portfolio A (for a total of H^). Thus H satisfies the 
recursion 

H k = (H k „ x - yV^ x )R + -yV^A + e k . (60) 
Error Propagation 

The replication error D k ~V k ~H k satisfies 

D k = RD k + e k (6i) 
where e k is uncorrelated with the market and 

e k = V k - V k - cov(V k , A)(A - 3)/oi (62) 

It follows that 

var(£ fc ) = var(T4) - cov(V kl A) 2 /o* A . (63) 
In terms of the lattice parameters this becomes 

~M = {p.d - P.) - ( ^XJ - Vtf. (64) 

If the dynamics of the variables A and B are stationary, the expression in brackets is 
constant, and hence the error variance is simply a constant times the square of the 
difference in the two successor V k &. 

LetU k ^lT- k D k . Then^ = ^ 1 + iF" fc e & . If S k is the variance of U T as seen 
starting at k } then 

«U = E[S k ] + iP^-*)var(6 Jb ) (65) 

with terminal condition S T = 0. This can now be. evaluated by a backward recursion 
in the lattice. 
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Method Steps 

Following is a summary of the steps of the lattice method for the discrete-time case. 
These steps are depicted in PIG. 7 and parallel those for the continuous-time case 
illustrated in FIG. 2. 

1. Step 700 is the set up. Given a random payoff G that depends on a non-traded 
variable B that evolves randomly, formulate a binomial lattice model. 

2. Step 710 determines an appropriate market representative A as either (a) a 
market asset most correlated to £ or (b) the Markowitz or market portfolio. 

3. Step 720 determines the risk-neutral probabilities by using the formula (56). 

4. Step 730 solves for the values of V on the lattice corresponding to the variable B 
using the recursion (57). Initiate with the terminal boundary condition V = G. 
The result at time 0 is the price. 

5. Optionally, step 740 determines the optimal replicating portfolio from H 0 - V 0 , 
and thereafter at step k-1 investing H h - X - in the risk free asset and 
7Vi_i in the asset most correlated with B. 

6. Optionally, step 700 determines the error variance by the backward recursion 
(65). 

IMPLEMENTATION 

The present invention may be implemented in various different ways on a computer. 
For the purposes of the present description, the term 'computer 1 is defined to include 
any electronic digital information processor, including hand-held calculators, personal 
data assistants, pocket personal computers, laptop computers, desktop computers, 
and so on. The computer may be programmed by providing a computer readable 
digital storage medium containing instructions to execute various methods of the 
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invention. The program or programs are then executed on the computer. Information 
associated with the variable x m may be provided to the computer manually by a user, 
received from a data network, or retrieved from a storage medium into computer 
memory. The program then retrieves the information associated with the variable 
x m from the computer memory and performs calculations in accordance with various 
methods of the invention. The results of the calculations may be used to display or 
otherwise communicate information to the user, or they may be stored or transmitted 
digitally for further processing or later display. 

More specifically, the implementation of methods of the invention may take various 
different forms. Following are representative examples of specific implementations. 

1. Numerical Solution. The straightforward way to implement the continuous-time 
method once the model is denned is by solution of the extended Black-Scholes 
equation (15). This can be carried out by standard numerical procedures, fol- 
lowing the methods developed for the ordinary Black-Scholes equation. The 
simple, one-dimensional version shown explicitly in (15) can be solved by a finite 
grid method: a large two-dimensional grid of points is denned with coordinates 
corresponding to x c and t and grid point values corresponding to the associated 
V values. In the simplest method, first-order derivatives of V are formed as 
(normalized) differences in grid point values. The second-order derivative is the 
second-order difference of grid point values. Boundary conditions are typically 
imposed at t = T, and at lower and upper values of x e . In practice, if the 
problem is defined for all values of x e with 0 < x c < co, an appropriate value 
for V(x e> t) at a high upper boundary curve x e (t) is assigned as a boundary 
condition in addition to the tenninal condition and the condition at x c = 0. It 
is possible to carry out the complete solution with a spreadsheet program such 
as Excel rurming on a general purpose computer. Experience has shown that 
the primary technical concern for simple methods such as this is that the time 
step between successive grid points be small in order to assure convergence. 
Small grid point widths in both x and t are used to get good accuracy. 
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For higher-order equations associated with extensions of the basic method, it is 
convenient to use a professional software package or one written in a program- 
ming language. The value of a call option on a non-tradable variable, but with 
a stochastically varying interest rate has been solved by the inventor using the 
Excel Visual Basic package. 

The same techniques may be applied to solve for the variance according to the 
partial differential equation (45). In simple experiments, the same grid size as 
used for finding V worked for rinding the non-hedgable variance. 

2. Explicit Solution. The explicit solution for call options (or for puts using the 
put-call parity formula) is easily solved by using either a table, or more likely, 
an approximate formula for values of the normal distribution. Again this can 
be carried out with Excel. It could be easily incorporated into a hand-held 
calculator. 

3. Discrete-time Version. The discrete time recursion is easily carried out by direct 
recursion. This method is relatively free from convergence issues, since both the 
model and the solution method are carried out with the same time steps (unlike 
the continuous-time case which is computed with a discrete process). 

4. Lattice Methods. The lattice method is easily implemented. In simple cases it 
may be carried out with a spreadsheet program running on a computer. The 
method becomes more challenging to implement when the model parameters 
are not constant, and when higher-order lattices and trees are used. 

5. Risk-Neutral Computation. Beginning with a discrete model, the risk-neutral 
version of the processes can be implemented computationally by computing the 
risk-neutral expected value of the payoff function. This can be carried out with 
backward recursion, simulation, or by constructing a lattice that has the risk- 
neutral probabilities on its arcs and using backward recursion on that lattice. 

6. Simulation. Simulation is a powerful method for implementing methods of this 
invention. Standard simulation packages can be used, such as a spreadsheet 
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program like Excel, modest simulation packages such as Crystal Ball, or more 
advanced statistical packages. At a professional level, special software that 
accounts for the financial structure as well as advanced simulation concepts 
would be used. 

7. Optimization. In many situations, on-going decisions can be made that influence 
the value of the payoff. For example, an American-style option allows the owner 
to exercise the option at any time before expiration, and selection of this exercise 
policy is an important component of the analysis of the option. Similarly, in a 
business venture there are opportunities to expand, contract, delay, and so forth. 
These policies can be found with the methods of this document in conjunction 
with standard methods for policy optimization employed for ordinary derivative 
theory. 
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CLAIMS 

What is claimed is: 

1. A computer-implemented method for pricing a financial derivative of a non- 
marketed variable x c , the method comprising: 

a) determining a market representative x m useful in detennining a value of the 
financial derivative; 

b) retrieving information associated with the non-marketed variable x e and the 
market representative x m ; 

c) calculating a solution to an equation involving a price of the financial derivative 
V(x ey t) defined as a function of x e and time t y wherein the equation comprises a 
coefficient involving the information associated with x e and x m ; 

and 

d) generating an output including the calculated price of the financial derivative. 

2. The method of claim 1 wherein the information associated with x e and x m 
comprises a drift rate of the non-marketed variable x e> and a drift rate of the market 
representative x m . 

3. The method of claim 1 wherein the information associated with x e and x m 
comprises variances of the non-marketed variable x c and the market representative x mt 
and a covariance between the non-marketed variable x e and the market representative 

4. The method of claim 1 wherein the coefficient involving the information as- 
sociated with x e and x m has the form ^ - /Mm*> - r), where ^ is a drift rate of 
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the non-marketed variable x ej is a drift rate of the market representative x m% r 
is an interest rate, and fl^ is a factor derived from a variance of the market repre- 
sentative x m and a covariance between the non-marketed variable x c and the market 
representative x m . 

5. The method of claim 1 wherein the equation is a modified Black-Scholes equa- 
tion. 

6. The method of claim 5 wherein the modified Black-Scholes equation is obtained 
from a standard Black-Scholes equation by replacing, in a term involving a first-order 
partial derivative of V(x et t) with respect to x e , a coefficient r, representing an interest 
rate, by a coefficient involving the information associated with x e and x m . 

7. The method of claim 1 wherein the equation is a discrete-time equation involv- 
ing V(x ei t) defined as a function of x e and discrete time points t = k. 

8. The method of claim 1 wherein the market representative x m comprises a 
marketed asset or combination of such assets that is approximately most correlated 
with the non-marketed variable x e . 

9. The method of claim 1 wherein the market representative x m comprises a 
combination of multiple marketed assets associated with market sectors most closely 
associated with the non-marketed variable x e . 

10. The method of claim 1 wherein the market representative x m comprises a 
marketed asset or combination of such assets that is approximately equal to an overall 
market portfolio. 
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11. The method of claim 1 further comprising calculating an optimal hedge. 

12. The method of claim 1 further comprising calculating a minimum variance of 
the error between an optimal hedge and the calculated price of the financial derivative. 

13. The method of claim 1 wherein the equation represents a risk-neutral dis- 
counted expected value of cash flows of the financial derivative. 

14. The method of claim 13 wherein a cash flow of the financial derivative is 
path-dependent. 

15. The method of claim 1 applied to derivatives of a set of non-marketed variables 
wherein the market representative x m comprises a combination of multiple marketed 
assets, each most-correlated with a different non-marketed variable in the set of non- 
marketed variables. 

16. The method of claim 1 wherein the calculated price of the financial derivative 
includes cash flows at an intermediate time and a terminal time. 

17. The method of claim 1 wherein drift rates, an interest rate, variances, and 
covariances of x t and x m either vary with time or are governed by stochastic processes. 

18. The method of claim 1 wherein the cash flow depends on marketed variables 
as well as non-marketed variables. 

19. The method of claim 1 wherein the equation involves additional non-marketed 
variables. 
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20. The method of claim 1 wherein the market representative is derived from a 
combination of multiple marketed variables, and wherein x e and the multiple mar- 
keted variables are governed by either geometric Brownian motion or alternative 
processes. 

21. A computer-implemented method of pricing a financial derivative of a non- 
marketed finite-state variable B y the method comprising: 

a) determining a market representative A associated with the non-marketed finite- 
state variable B\ 

b) calculating risk-neutral probabilities for the non-marketed finite-state variable 
B using a binomial lattice model associated with the non-marketed finite-state vari- 
able B and the market representative A; 

c) calculating values of a price function V defined on the lattice corresponding to 
the variable B\ 

and 

d) generating from the calculated values of the price function V an output includ- 
ing a calculated price of the financial derivative. 

22. The method of claim 21 wherein the market representative A is determined to 
be approximately equal to at least one of a Markowitz portfolio, a market portfolio, 
and a market asset most correlated to the non-marketed finite-state variable B. 

23. The method of claim 21 further comprising calculating an optimal hedge. 

24. The method of claim 21 further comprising calculating a Tr^mum variance of 
the error between an optimal hedge and the calculated price of the financial derivative. 
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26. The method of claim 21 wherein a cash flow of the financial derivative is 
path-dependent. 

26. The method of claim 21 wherein the binomial lattice model comprises time- 
dependent lattice parameters of the variables A and B. 
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